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Abstract. We study concentrated colloidal suspensions, a model system which has a glass transition. The non-equilibrium 
nature of the glassy state is most clearly highlighted by aging - the dependence of the system's properties on the time elapsed 
since vitrification. Fast laser scanning confocal microscopy allows us to image a colloidal glass and track the particles in three 
dimensions. We analyze the static structure in terms of tetrahedral packing. We find that while the aging of the suspension 
clearly affects its dynamics, none of the geometrical quantities associated with tetrahedra change with age. 



INTRODUCTION 

If a liquid is rapidly cooled faster than it can crystallize, 
it can form an amorphous solid also known as a glass. 
This occurs below a temperature known as the glass 
transition temperature, To. As To is approached from 
above, the viscosity of the liquid increases by many 
orders of magnitude over a small range in temperature. 
However, a diverging structural length scale that would 
also grow by orders of magnitude and thus explain this 
thickening has not been identified [1-5]. 

A related phenomenon is that of aging [6, 7], the ob- 
servation that a glassy sample is not in equilibrium and 
thus its properties depend on the age of the sample. Were 
the system able to equilibrate at or below T g it would 
assume the lowest energy configuration: a crystal. How- 
ever upon a temperature quench to T < T g the system's 
dynamics are severely slowed and its structure is frozen 
in an amorphous configuration. Slight thermal motions 
do survive this temperature quench and slowly allow the 
system to change configurations albeit in a non-ergodic 
fashion. The time scale for rearrangements grows with 
the age of the sample; the older the sample, the longer 
the time needed for a further change. Similarly to the 
problem of the glass transition, a changing time scale 
for aging dynamics is suggestive of an underlying chang- 
ing structural length scale. One conjecture is that as the 
glass ages it optimizes its packing. Specifically the length 
scale over which some "optimal packing" is achieved 
might grow during aging, thus slowing the dynamics of 
the sample, which in turn would slow the aging process 
itself. 

Suspensions of colloids have been studied as a model 
glass-forming system with great success. Colloidal sus- 
pensions consist of solid particles in a liquid, and exhibit 
a glass transition as their concentration is increased. In 



such systems, the control parameter is not temperature 
but packing fraction, 0, and a glass is obtained once <p 
increases above 2 0.58 [8-11]. Previous work with 
colloids [12, 13] has shed some light on the aging pro- 
cess showing, for example, that aging is both spatially 
and temporally heterogeneous. However a detailed un- 
derstanding of the microscopic mechanisms involved is 
still lacking. 

In this work we analyze aging in a colloidal glass from 
a geometrical point of view. In particular we focus on the 
idea of geometrical frustration [14, 15]. In hard sphere 
systems the free energy is exclusively determined by en- 
tropy and, in order to minimize the former at a given 
temperature, the system has no choice but to maximize 
the latter. This is evident already at a packing fraction 
of ^crystallization ~ -494 where the system crystallizes, de- 
spite the lack of inter-particle interactions, sacrificing 
some configurational entropy to gain vibrational entropy 
instead [16]. This is achieved by re-distributing the avail- 
able space homogeneously so that all particles have more 
local free volume and thus room to move about their av- 
erage lattice position. 

The most efficient packing of four spheres in three 
dimensions is to pack them into a regular tetrahedron. 
The effective volume fraction of four such spheres is 
^tetrahedron ~ 0.78. However regular tetrahedra don't tile 
3-D space, and thus the most efficient space filling con- 
figuration is a hexagonally closed packed crystal with 
ss 0.74. So while tetrahedral packing could locally 
maximize vibrational entropy, the constraint that the 
sample fill space induces crystallization. This frustration 
between local and global packing optimization has been 
invoked as a possible cause for the glass transition in sim- 
ple liquids [14, 15]. 

We use fast laser scanning confocal microscopy to ob- 
serve micron-sized particles in a colloidal glass in three 
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dimensions. Specifically we study the static structure of 
the glass in terms of tetrahedral packing and analyze the 
evolution of some geometrical properties as the sample 
ages. While the dynamics change dramatically over the 
course of our experiment, we find no change in the struc- 
tural properties, as characterized by tetrahedra arrange- 
ments. 



EXPERIMENTAL DETAILS 

We use poly(methyl methacrylate) (PMMA) particles of 
diameter d = 236jim and polydispersity ~ 5%. The par- 
ticles are sterically stabilized by a thin layer of poly- 12- 
hydroxystearic acid, dyed with rhodamine 6G and sus- 
pended in a mixture of cyclohexylbromide and decalin 
(85:15 by weight) that closely matches both the density 
and index of refraction of the particles. This is done to 
minimize sedimentation effects and to improve imaging. 
The solvent's viscosity at 295 K is 77 = 2.25 mPa-s which 
sets the time required to diffuse a distance d in the dilute 
limit to gj = Tdiff - 23.5 where D = 

We use a confocal microscope to acquire three- 
dimensional images of a viewing volume of 62 fim x 
58 /ira x 10 jJ,m at a rate of 1 frame every 26 s. Note that 
the sample is glassy - at these densities the colloids move 
very slowly. The particles can be tracked even though the 
frame rate is comparable to Tdiff. The viewing volume 
typically contains <~ 2500 particles. We focus at least 60 
jUw away from the cover slip of the sample chamber to 
avoid wall effects [19, 20]. We identify particles with a 
horizontal accuracy of 0.03 [im and a vertical accuracy 
of 0.05 jJ,m, and track them in three dimensions over the 
course of the experiment [17, 18]. 

The control parameter for colloidal phase behavior 
is the sample volume fraction 0. While the rhodamine 
imparts a slight charge upon the particles, their phase 
behavior, 0f ree ze = 0.38 and me i t = 0.42, is similar to 
that of hard spheres (0f ree ze = 0.494 and me i t = 0.545). 
We observe a glass transition at g ss 0.58, in agreement 
with what is seen for hard spheres. We examine samples 
with volume fractions in the range <p s=s 0.58 — 0.62 and 
find qualitatively identical results. These samples form 
small crystals which nucleate at the cover slip, but do not 
form crystals within the bulk of the sample even after 
several weeks. Here we report results based on a sample 
at w 0.62. 

The sample is initialized by vigorous stirring. The sub- 
sequent particle dynamics are reproducible after this stir- 
ring, and depend only on the waiting time f w since the 
start of the aging process. Within a minute of ending the 
stirring, transient flows within the sample greatly dimin- 
ish, and the particles move slowly enough to be iden- 
tified and tracked. This defines our initial time f w = 0, 
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FIGURE 1. Aging mean squared displacement for a col- 
loidal glass at (j) ~ 0.62. The three curves represent three dif- 
ferent ages of the sample. A : t w = OTjjff, x : f w = 50Tdiff and 
O : f w = ISOTjiff. The dashed line has a slope of 1 and repre- 
sents diffusive behavior, not seen in this glassy sample. 

or age zero, for the sample. The results below are not 
sensitive to variations of this choice. During the experi- 
ments, no crystallization is observed within the viewing 
volume. Note that in many aging studies, the initial sam- 
ple is prepared by a temperature quench (corresponding 
to a rapid increase of the volume fraction in our exper- 
iments). The initial conditions in our experiments cor- 
respond to a shear-melted sample; the volume fraction 
remains constant. 

Henceforth all numerical values will be expressed in 
terms of particle diameter d and characteristic diffusion 
time T d iff. 

RESULTS 

In order to determine whether our sample is aging we 
take a long (<~ 350Tdiff) data set and split it in three 
time windows as follows: [0 — 50Tdiff], [50 — 150Tdiff] and 
[150— 350Tdiff]. Since the sample is not disturbed dur- 
ing observation this is equivalent to running three exper- 
iments with f w = 0, 50, and 150Tdiff from the sample ini- 
tialization. We then measure the mean squared displace- 
ment for the three data sets averaged over all particles 
and over all initial times within a given window. As can 
be seen in Fig. 1 the data clearly show a slowing down of 
the dynamics. For example, as the sample ages (f w goes 
from to 150Tdiff) the time particles take to diffuse, say, a 
distance of 0.02c/ 2 goes from approximately At s=a 30Tdiff 
to At w 90T d iff. 

While the dynamics are clearly changing with the age 
of the sample, we wish to consider the static structure as 
well. We start by calculating the pair correlation function 
g(r). The resulting curve is plotted in Fig. 2, which has 
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FIGURE 2. Pair correlation function for a sample with <j> ~ 
0.62. The dashed segment indicates the range of interparticle 
distances used as a definition of nearest neighbor. 



been calculated from data averaged over all ages of the 
sample; no change in g(r) is seen as the sample ages. The 
first peak of g(r) is at r = 1 .04d, slightly displaced from 
the ideal hard-sphere position of r = d. This is mainly 
due to a slight charge induced by the dyeing process, 
and may also be due to some uncertainty in the exact 
value of the average particle size d which we estimate 
to be 2%. Another effect of this soft potential is that 
g(r) for r < d does not fall to zero as sharply as it 
would for perfect hard spheres. Figure 2 also shows the 
characteristic double peak around r = 2d seen in many 
glassy systems. 

To further characterize the geometry of the colloidal 
glass we identify tetrahedra formed by the particles as 
follows. Two colloids are labeled as nearest neighbors if 
their distance is within the first peak of g(r), namely, if 
0.74c/ > r > 1 .38of as shown in Fig. 2. The upper bound 
corresponds to the first minimum of g(r) while the lower 
bound is set to eliminate some nearest neighbor pairs that 
are excessively close to each other as a result of occa- 
sional errors in particle tracking. The results presented 
here are insensitive to variations in the limits for r and to 
the use of Delaunay triangulation as an alternative near- 
est neighbor finding method. 

A tetrahedron is then defined as a quadruplet of mutu- 
ally nearest neighbor colloids. The six inter-particle dis- 
tances, or "bond" lengths, b\ with i — 0...5 are measured. 
We then assign to each tetrahedron a "looseness" as mea- 
sured by the average bond length b and an "irregularity" 
as measured by the standard deviation of the six bond 
lengths (7b and defined by = Lf=o(^ ~bi) 2 . Note that 
the looseness has the same limits as the nearest neigh- 
bor distance (b £ [0.74d, 1.38c/]) while the irregularity is 
bound by zero (for an ideal tetrahedron) and by the max- 
imum of the standard deviation of six numbers drawn 
from the interval [0.74, 1.38], namely rj b & [0,0.35c/]. 





FIGURE 3. Computer generated picture of a very regular 
tetrahedron (Ob = 0.02d, left) and a very irregular one (ob = 
0.22c/, right). The spheres are drawn to scale (d = 2.36/im). 



Figure 3 shows a computer rendered image of a very 
regular and a very irregular tetrahedron with Oi, = 0.02a! 
and 0\, = 0.22c/ respectively. Clearly our definition of 
what constitutes a tetrahedron combined with the range 
we set for nearest neighbor distances allows for fairly 
irregular, quasi planar tetrahedra; though setting r max = 
1.38c/ < \/2d avoids the possibility of a perfectly planar 
configuration. 

Since aging is the evolution of the system's proper- 
ties over times much longer than Tdiff we analyze these 
two static quantities, b and Ob, as the sample ages. If ag- 
ing is a result of thermal fluctuations slowly trying to re- 
solve the frustration we referred to above, then one might 
expect the geometrical characteristics of the tetrahedra 
forming the sample to evolve with time. We therefore 
tracked the values of (b) and (rjb) averaged over all tetra- 
hedra in each frame as the sample ages, and show this in 
Fig. 4. These two average quantities remain essentially 
constant over the duration of the experiment, (b) varies 
by less than 0.5%, reflecting the fact that the packing 
fraction <j) remains constant. One could envisage that a 
change in the value of (ob) would allow the sample to 
find a configuration with lower energy. For example, a 
"better packed" configuration might have fewer irregular 
tetrahedra, as is the case for a crystal. Such a configura- 
tion might allow each particle to have more local room to 
move, thus increasing its vibrational entropy. This how- 
ever is not the case, as is seen in Fig. 4, which shows that 
(Ob) decreases by less than 0.5% throughout the experi- 
ment. While a decrease is consistent with our argument, 
this decrease is not reproducible in other samples, and in 
fact we find in all cases that (<7b) varies by at most 1% 
with no preference for increasing, decreasing, or merely 
fluctuating. 

As neither of these average static quantities evolve as 
the sample ages, we consider a more complete picture by 
calculating the probability distribution functions P(b/d) 
and P(&b/d). For each value of the f w the distributions 
were calculated by averaging over all tetrahedra and all 
times in each time window. 

The distribution of average bond lengths P(b/d) is 
shown in Fig. 5 and is peaked at b/d « (b/d) = 1.09. 
P(b/d) is approximately a Gaussian with a width of 0.04, 
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FIGURE 4. Plot of the average tetrahedral "loosness" (b) 
(top) and average tetrahedral irregularity (oi,) (bottom) versus 
time. The average is taken over all tetrahedra at each time. 



although the distribution is slightly asymmetric with a 
bias toward larger values of b. Thus compact tetrahe- 
dra (b/d < 1.05) are less prevalent, along with loose 
tetrahedra (b/d > 1.13). More significantly, the distri- 
bution P(b/d) is unchanged during aging, as is seen by 
the excellent agreement between the distributions plot- 
ted with the different symbols. Thus, while both compact 
and loose tetrahedra are less likely, their relative preva- 
lence remains unchanged as the sample ages. This seems 
somewhat surprising, as entropy might favor a more stan- 
dardized tetrahedron size. Compact tetrahedra have less 
room for their particles to move; the opposite is true for 
loose tetrahedra. The system could maximize its overall 
vibrational entropy by sharing free volume more evenly. 
However, this is not reflected in the data, and may in- 
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FIGURE 5. The distribution function P(b/d) for tetrahedral 
bond lengths ("looseness") for three ages. A : f w = 0, x : f w = 
50T(jiff, and O : f w = lSOT^ff . The solid line is a reference curve 
computed by considering random bonds based on g(r), without 
the constraint that such bonds form a tetrahedron. 
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FIGURE 6. The distribution function P(o b /d) for tetrahe- 
dral "irregularity" for three ages: A : f w = 0, x : ? w = 501^, 
and O : t w = 150T(jiff . The solid line is a reference curve com- 
puted by considering random bonds based ong(r), without the 
constraint that such bonds form a tetrahedron. 



dicate that random packing requires both compact and 
loose tetrahedra in order to tile space. Thus, the age of 
the sample is not detected by P(b/d). 

Now consider the distribution of tetrahedral irregular- 
ities P(a b /d) plotted in Fig. 6. The most striking feature 
is a broad maximum around 0\,/b w 0.1. This is con- 
gruent with the amorphous quality of a glassy packing 
where, unlike in a crystal, one might expect to find a great 
variety of tetrahedra, from quite regular to quite irregu- 
lar ones. However, it is also interesting to note that the 
probability of finding a tetrahedron with an irregularity 
greater than 0.2 is effectively zero while the theoretical 
maximum for Ob/d is 0.35. Also, while not very regular, 
the shape of P{a b /d) is qualitatively reproducible with 
other samples above <j) g . 

The striking agreement in the tetrahedral properties 
between different ages of the sample suggests that there 
may be a simple underlying explanation for these distri- 
butions. One hypothesis is that the shape of these distri- 
butions is dictated by the statistics of choosing six bond 
lengths bi at random from the overall distribution P(b{) 
in an amorphous configuration. P(bi) is closely related to 
g(r) by P(r) = r 2 g(r). By randomly picking many sets 
of six bi, without the constraint that they form a tetra- 
hedron, we can compute reference distributions P'(b/d) 
and P'(o\,/d). These are shown as solid lines in Figs. 5 
and 6. While these are similar to the experimental data, 
deviations are clearly seen. In Fig. 5, the solid curve 
shows a slight bias toward larger values of b/d as com- 
pared with the experimental data. In Fig. 6, the peak 
of the distributions differs noticeably between the ex- 
perimental values and the reference curve. While these 
differences are not extreme, they are outside of the un- 
certainty of the data, and thus indicate that tetrahedral 
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properties do possess additional information beyond just 
nearest neighbor bond lengths [which in turn are related 
tog(r)]. 

CONCLUSION 

We have examined an aging colloidal glass by looking 
at static structural characteristics of tetrahedral packing. 
We find that there is no correlation between the age 
of a glass and its static structure as evaluated by these 
quantities. Nevertheless, our sample's "awareness" of 
its age must be encoded in the particle positions, and 
thus in some feature of the packing of the particles. 
We have checked other geometrical quantities based on 
tetrahedral structure, such as volume or surface area, and 
none of these show aging either. Tetrahedral properties 
appear not to reflect the age and thus a static measure of 
a sample's age remains elusive. 
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